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Abstract 

Within the description of stochastic differential equations it is argued that the existence of 
Boltzmann-Gibbs type distribution in economy is independent of the time reversal symmetry in 
econodynamics. Both power law and exponential distributions can be accommodated by it. The 
demonstration is based on a mathematical structure discovered during a study in gene regulatory 
network dynamics. Further possible analogy between equilibrium economy and thermodynamics 
is explored. 

PACS numbers: 

89.65.Gh Economics; econophysics, financial markets, business and management; 
05.10.Gg Stochastic analysis methods (Fokker-Planck, Langevin, etc) 
05.20.-y Classical statistical mechanics. 
87.23.Ge Dynamics of social systems 
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I. INTRODUCTION 



Carefully and extensive analysis of real economic and financial data have revealed vari- 
ous exponential and power-law distributions regarding to money, income, wealth, and other 
economic quantities in ancient and modern social societies . A remarkable anal- 

ogy between the economic system and a thermodynamical system has been revealed in a 
recent study [4]. Using a detailed microdynamical model with time reversal symmetry, it 
was demonstrated that a Boltzmann-Gibbs type distribution exists in economic systems. 
Indeed, ample empirical analysis support such suggestion jj, O]. Nevertheless, different mi- 
crodynamical models lead to apparently different distributions 0, Q]. Those distributions 
are supported by empirical data, too. The nature of such difference may reveal the differ- 
ence in corresponding economic structure. For example, this difference has been tentatively 
attributed to the role played by time reversal symmetry in microdynamical models 0]. 

Following the tradition of synthesizing interdisciplinary knowledge, such as from biology, 
physics, and finance p], in this letter we argue that irrespective of the time reversal symmetry 
the Boltzmann-Gibbs type distribution always exists. A broader theoretical base is thus 
provided. The demonstration is performed within the framework of stochastic differential 
equations and is based on a novel mathematical structure discovered during a recent study in 
gene regulatory network dynamics |7[. In the light of thermodynamical description, possible 
explanations for the origin of the difference in various empirical distributions are proposed. 

II. BOLTZMANN-GIBBS DISTRIBUTION IN FINANCE 

Stochastic differential equations or Langevin equations and their corresponding Fokker- 
Planck equations have been shown to be a useful modelling tool in economy 03,0, EJ- One 
of the best examples is the Black-Scholes formula in option pricing theory [8|. This kind 
of mathematical formulation provides a direct connections between the microdynamics and 
the station^ state and has been nsed to g ene r ate various dl8 t ribu t io n >aws !MQg. 

Specifically, the stochastic differential equation may take the following form: 

q = f(q)+iV 7 (q)£(t), (1) 

where f and q are n-dimensional vectors and f a nonlinear function of the state variable q. 
The state variable is the quantity to specify the economic system. It may be the money, 



the income, or other suitable indices. The noise £ is a standard Gaussian white noise with / 
independent components: (&) = 0, (£i(i)£j(t')) = 26ij6(t — t'), and i, j = 1,2, In Eq.(l) 
we explicitly factorize out the pure noise source and the state variable dependent part for 
the convenient of later description. 

The specification of the noise in Eq.Q is through the n x n diffusion matrix -D(q) by the 
following matrix equation 

Ni{ci)NJ (q) = e D(q) , (2) 

where Nj is an n x I matrix, Nj is its the transpose, and e a nonnegative numerical constant 
to keep tract of the noise. It plays the role of temperature in thermodynamics. According to 
Eq.(2) the n x n diffusion matrix D is both symmetric and nonnegative. For the dynamics 
of state variable q, all what needed from the noise is the diffusion matrix D. Hence, it is 
not necessary to require the dimension of the noise vector £ be the same as that of the state 
vector q and to require more specific knowledge of n x I matrices {Nj} beyond Eq.(2). 

It is known that even in situations that Eq. (Q) is not an exact description, and perhaps it 
would never be in a rigorous sense in economy, it may still serve as the first approximation 
for further modelling [l J2 L Indeed, it has been empirically verified to be a rather accurate 
description in economy a, |9|, [lCJ, [ll ] . Because the energy function or Hamiltonian has played 
a dominant role in equilibrium physics processes, the crucial question is whether or not a 
similar quantity exists in a more general setting. In the following we present an argument 
leading to the positive answer. 



There exists several ways to deal with the stochastic equations equation in the form 



of 



q.(l) and (2). The most commonly used are those of Ito and Stratonovich methods 

as 

, |12J. However, with those methods the connection between the existence of energy 
function like quantity in Eq.(l) and the stationary distribution is not clear when the time 



reversal symmetry is broken [12j. The difficulty for finding such potential function can be 
illustrate d by the fact that usually _D _1 (q)f (q) cannot be written as the gradient of a scalar 
function 12| in the absence of detailed balance condition or in the broken time reversal 
symmetry. This will become precise as we proceed. 

During a recent study of the robustness of the genetic switch in a living organism it 
was discovered that Eq.(l) can be transformed into the following form, 

[A(q) + C(q)]q = «9 q 0(q) + N n {q)£{t) , (3) 
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where the noise £(£) is from the same source as that in Eq.(l). Here we tentatively name 
A(q) the adaptation matrix, C(q) the conservation matrix, and the scalar function <f)(q) the 
fortune function. The gradient operation in state space is denoted by <9 q . The adaptation 
matrix A(q) is defined through the following matrix equation 

N II (q)NJ I (q) = eA(q) 1 (4) 

which guarantees that A is both symmetric and nonnegative. The nxn conservation matrix 
C is antisymmetric. We define 

A(q) + C(q) = l/[D(q) + Q(q)] = M(q) . 

with the nxn matrix M is the solution of following two matrix equations Q|. The first 
equation is the potential condition 

q x [M(q)f (q)] = , (5) 

which gives n(n — l)/2 conditions [the wedge product for two arbitrary vectors Vi and v 2 
in n-dimension: [vi x v 2 ]ij = v u V2j — v\jV<n , i, j = 1,2, ...,n ]. The second equation is the 
generalized Einstein relation between the adaptation and diffusion matrices in the presence 
of conservation matrix 

M(q)D (q) M T (q) = ±[Af(q) + ikT(q)] . (6) 

which gives n(n + l)/2 conditions The fortune function 0(q) is connected to the 

deterministic force f(q) by 

9 q 0(q) = M(q)f(q) . 

7or simplicity we will assume det(A) ^ in the rest of the letter. Hence det(M) ^ 
1^ . Thus, the adaptation matrix A, the conservation matrix Q and the fortune function 
<p in Eq.(3) and (4) can be completely determined by Eq.(l) and (2). The breakdown of 
detailed balance condition or the time reversal symmetry is represented by the finiteness of 
the conservation matrix 

C(q) + , (7) 

or equivalently Q ^ 0. The usefulness of the formulation of Eq.(3) and (4) is already man- 
ifested in the successful solution of outstanding stable puzzle in gene regulatory dynamics 



and in solving two fundamental controversies in population genetics 
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A few remarks on Eq.(3) are in order. In the light of classical mechanics in physics, 
Eq.(3) is in precisely the form of Langevin equation. The fortune function <ft corresponds 
to the potential function but opposite in sign to reflect the fact that in economy there is a 
tendency to seek the peak or maximum of fortune. The adaptive matrix A plays the role 
of friction. It represents adaptive dynamics and is the dynamical mechanism to seek the 
nearby fortune peak. The conservation matrix C plays the role analogous to a magnetic 
field. Its dynamics is similar to that of the Lorentz force, hence conserves the fortune. As 
in classical mechanics, the finiteness of the conservation matrix C breaks the time reversal 
symmetry. 

It was heuristically argued [k| and rigorous demonstrated [16] that the stationary dis- 
tribution p(q) in the state space is, if exists, 

Therefore, the fortune function acquires both the dynamical meaning through Eq.© 
and the steady state meaning through Eq.({8|). Specifically, in the so-called zero-mass limit 
to differentiate from Ito and Stratonovich methods, the Fokker-Planck equation for the 
probability distribution p(q, t) takes the form [3| 

d t p(q, t) = d^M- 1 (q) [ed q - <9 q 0(q)]p(q, *) ■ (9) 

Here dt is a derivative with respect to time and <9 q represents the gradient operation in state 
space. We note that Eq.fJHJ) is a stationary solution to Eq.(jOJ) even it may not be normalizable, 
that is, even when the partition function Z = J d n q p(q) is ill-defined. Again, we emphasize 
that no time reversal symmetry is assumed in reaching this result. This completes our 
demonstration on the existence of the Boltzmann-Gibbs distribution in economy. 

Using M _1 (q) = D(q) + O(q) and f(q) = [D(q) + Q(q)]<9 q 0(q), Eq.© can be rewritten 
in a more suggestive form [16] 

d t p(q 7 t) =^[eD(q)9 q + e(^Q(q))-f(q)]p(q,t) . (10) 

It is clear that in the presence of time reversal symmetry, i.e. Q = 0, one can directly read 
the fortune function from above form of Fokker-Planck equation as <9 q 0(q) = D~ 1 (q)i (q). 
For the sake of completeness, we list the Fokker-Placnk equations corresponding to Ito 
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and Stratonovich treatments of Eq.(l) and (2) 



121: 



d t pi(q,t) =J2 d <i 



j'=i 



pi(q,t) , (Ito) 



and 



n £ 

■/((ql + EE^WWi) 

j=i ^=1 



p s (q,t). (5*r) 



(11) 



;i2) 



The connection of fortune function with both dynamical (Eq.(l) and (2)) and stationary 
state is indeed not clear in above two equations. Nevertheless, it has been shown [16] that 
there are corresponding fortune functions, adaptation and conservation matrices to Eq.([ll|) 
and (|12j). We point out here that when the matrix Ni is independent of state variable, 
Eq. ffTTj) and (fT2"|) are the same but may still differ from Eq.(|5|). because the gradient of the 
antisymmetric matrix Q(q) may not be zero. This last property shows that the time reversal 
symmetry is indeed important. 



III. TWO EXAMPLES 

The Fokker-Planck equation used by Silva and Yakovenko has the form: 

d t psy(q, t) = d q [a(q) + d q b(q)]p sy (q, t) , (13) 

and Fokker-Planck equation used by Bouchaud and Mezard [3] has the form 

dtpbmiq, t) = d q [(J(q - 1) + a 2 q + a 2 qd q q}p bm (q, t) . (14) 

They are all in one dimension. We immediately conclude that the conservation matrix C, 
equivalently Q, is zero, because there is no conservation matrix in one dimension. In accor- 
dance with the definition in the present letter, which is consistent with that in nonequilibrium 
processes jl2| . the dynamics described by above two equations can be effectively classified 
as time reversal symmetric. 

Rewriting them in symmetric form with respect to the derivative of state variable q as in 
Eq. ffTUj) . we have 

d t psy{q, t) = d g [a(q) + {d q b{q)) + b(q)d q ]p sy {q, t) , (15) 

and 

d t pbm(q, t) = d q [(J(q - 1) + o 2 q + a 2 q + cr 2 qd q ]p bm (q, t) . (16) 



The corresponding wealth functions can be immediate read out as 



4>bm{q) 




11 ( j\ J i ( J\ , 

--- 2 + — lng + — - + 2 + — lng 

1 a \ a z J o On V <j ) 



(17) 



(18) 



They are exactly what found in Ref.[ja|] and [|3|]: the first one corresponds to an exponential 
distribution and the second one a power law distribution according to the Boltzmann-Gibbs 
distribution Eq.(jSJ). 

IV. ENSEMBLES AND STATE VARIABLES 

Having defined a precise meaning of time reversal symmetric and have demonstrated that 
the Boltzmann-Gibbs distribution even in the absence of time reversal symmetry, we explore 
further connection between econodynamics and statistical physics. 

There are two general types of situations which would generate different distributions in 
statistical physics and thermodynamics. The first one is to link to constraints on the system 
under various conditions. In statistical physics such constraints are described by various 
ensembles and free energies. For example, there are canonical and grand-canonical ensem- 
bles. There are Gibbs and Helmholtz free energies, entropy, enthalpy, etc. Those ensembles 
have their characteristic distributions. It would be interesting to know the corresponding 
situations in economy. 

Even with a given constraint, the form of distribution depends on the choice of state 
variable. For example, for ideal gas model, the distributions are different if views from the 
kinetic energy and from velocity. Hence, there is a question of appropriate state variable 
for a given situation, with which the physics becomes particular transparent. It would 
be interesting to know what be the appropriate variables to describe an economic system. 
Within this context, the difference between what discovered by Draeulecu and Yakovenko 
Q and Boachaad aad Mezarf S is perhaps mOI e due to the dinereace ia choice, of state 
variables, because it seems they are describing the same situation of same system under 
same constraints. 
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A remark on terminology is in order. It was demonstrated above that regardless of 
the time reversal symmetry the Boltzmann-Gibbs distribution, Eq.fjHJ), exists. The fortune 
function has an additional dynamical meaning defined in Eq.fjHJ). Both exponential and 
power law distribution can be represented by Eq.fjHJ). In fact, it is well known that power 
law distributions exist in statistical physics. A nontrivial example is the Kosterlitz-Thouless 



transition |17[. Thus it does not appear appropriate to call the power law distribution non- 
3oltzmann-Gibbs distribution. Such a terminology confusion was already noticed before 

fl n 

In the view of the dominant role of entropy in Kosterlitz-Thouless transition [121], the 
ubiquitous existence of power law distribution in economy may suggest that the entropy 
effect is rather important in econodynamics. This may corroborate with the suggestion of 
"superthermal" in economy 



V. CONCLUSIONS 

In this letter we demonstrate that the existence of Botlzmann-Gibbs distribution in fi- 
nance is independent of time reversal symmetry. Both power law and exponential distribu- 
tions are within its description. In analogous to similar situation in statistical physics, the 
differences among those distributions discovered empirically in economy are likely the result 
of different choices of state variables to describe the same system in econodynamics. 

This work was supported in part by USA NIH grant under HG002894. 
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